Epidemic spreading with immunization rate on complex networks 
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Abstract 

We investigate the spread of diseases, computer viruses or information on complex networks and also 
immunization strategies to prevent or control the spread. When an entire population cannot be immunized 
and the effect of immunization is not perfect, we need the targeted immunization with immunization rate. 
Under such a circumstance we calculate epidemic thresholds for the SIR and SIS epidemic models. It is shown 
that, in scale-free networks, the targeted immunization is effective only if the immunization rate is equal to 
one. We analyze here epidemic spreading on directed complex networks, but similar results are also valid for 
undirected ones. 



1. Introduction 

The spread of diseases, computer viruses or infor- 
mation in complex networks such as social networks, 
the WWW and Internet has attracted researchers' at- 
tention in various fields. We refer to [2, 3, 8] and ref- 
erences therein for abundant examples of complex net- 
works. Also preventing the spread has extensively been 
investigated. A typical method to prevent or control 
an epidemic spread is immunization such as vaccination 
[5]. Generally, it is difficult or impossible to perfectly 
immunize an entire population, so the object or target 
of immunization must be prescribed. However in turn, 
all of target individuals may not be immunized due to, 
for example, lack of plenty of time and resources, other- 
wise because some are overlooked. We describe this sit- 
uation as "imperfect immunization to targeted nodes". 
Nodes are sometimes called vertices or individuals. In 
another aspect, vaccination or a treatment cannot pro- 
vide sufficient immunity upon target nodes, which we 
call "insufficient immunity". Therefore, we have to in- 
troduce a parameter a (0 < a < 1) called immunization 
rate for such both cases. As in [12] we assume that a 
is a constant throughout. 

In [10] some immunization schemes are evaluated in 
the SIS model on undirected networks and [4] also of- 
fers other schemes for immunization, and [12] provides 
targeted immunization strategies in the SIR model. All 
those studies are devoted to undirected. In undirected 
networks links (or edges) do not have directions or they 
are bidirectional by nature. 

However, one can find many directed complex net- 
works [2, 3, 8] in nature, society, and artificial struc- 
tures; food webs, phone-call networks, the WWW, etc. 
The indegree of a node is the number of incoming links 
into the node and the outdegree is that of outgoing 
links emanating from it. If the indegree and outdegree 
distributions of a directed network follow power laws: 

P(k) oc k-' 1 and Q(£) oc 

for indegrees k and outdegrees £, and if the exponents 



satisfy 2 < 7,7' < 3, then the network is called scale- 
free. Undirected networks make no distinction between 
both degree distributions. 

The direction plays an important role in the study 
of epidemic spreading on a directed network. Diseases, 
computer viruses or information are transmitted to 
other nodes (or vertices) through outgoing links and 
a node is infected by incoming links. We derive the 
critical infection rate or the threshold, above which a 
disease spreads in the networks and below which it dies 
out. 

First we prescribe a target set of immunization in 
terms of indegrees and outdegrees on directed networks. 
We derive the thresholds for the SIR and SIS epidemic 
dynamics involving the target set and an immunization 
rate a. Actually they turn out the same for both. In 
the SIR model, the averaged fraction of nodes that are 
ever infected until the disease dies out is also given. 
Finally, based on thresholds, we will discuss whether 
immunization strategies are effective or not in scale- 
free networks. In particular, all targeted immunization 
schemes are shown to be effective only if a = 1, which 
was assumed in [10] from the beginning. 

Since undirected networks are simpler than directed 
ones, we immediately see that similar results hold on 
undirected networks by the arguments presented here, 
and they are also partially found in [4] for the SIS 
model and in [12] for the SIR model. 

2. Models of immunization in complex 
networks 

In this section we introduce preliminary definitions 
and model epidemic spreading with immunization rate. 
Our model is based on [11], where epidemic spread- 
ing with no immunization is discussed on directed net- 
works. First we treat the SIR model and later proceed 
to the SIS model which is regarded as a special case of 
the SIR model. 

In the SIR model, nodes of the network are di- 
vided into the following three groups regarding infec- 
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tion states ([5], [7, Chap. 10]): Susceptible (S), Infected 
(I) and Removed (R). In order to deal with immuniza- 
tion such as vaccination, we add one more state: Im- 
munized. Outgoing links emanating from immunized 
nodes are harmless. 

Hereafter we will denote a susceptible node by an 
S-node etc., for short. An S-node becomes infected at 
a rate of A (0 < A < 1). The parameter A is the infec- 
tion rate, for which we will derive the critical value for 
an outbreak of a disease, a computer virus, etc. The 
disease can be passed from I-nodes to S-nodes follow- 
ing only the direction of directed links. R-nodes have 
cither recovered from the disease or died and so they 
cannot pass the disease to others. An I-node becomes 
an R-node at a rate S (0 < S < 1) and, without loss of 
generality, we will set S = 1. 

Let p(k,£) be the joint probability distribution of 
nodes with indegree k and outdegree £, and let us de- 
note the marginal distributions by 

p(fc) = 5>(M), Q(e) = J2p( k > £ ) 

i k 

and the averages by 

k,t k 

(k£) =Y,k£p(k,£). 

k,e 

We denote by T the object or target set for immu- 
nization that is characterized in terms of indegrees and 
outdegrees, T being the complement of T that is not 
an object of immunization. The notation (k, £) e T 
means that the population of all nodes with indegree k 
and outdegree £ is an object of immunization. In case 
of T = {(k,£)\k > K}, for example, the population of 
all nodes with indegrees exceeding a number K is col- 
lectively an object of immunization. The summation 
restricted to all (k,£) in T will be denoted by and 
similarly for T by J2t- Furthermore, the average of the 
product k£ over T is denoted by (k£) T : 



(k£) 7 



J2 k£p(k,£) = Y, ke P( k > £ y 



(k,e)er 



In this paper we assume that the immunization rate 
a (0 < a < 1) is a constant as in [12], although it can 
be dependent on k and £ as well. In the following two 
models A and B, remark that a has different mean- 
ings, althought the same notation is used. Within the 
population of all nodes with indegree k and outdegree 
£, variables Sk,e(t), Pk,e(t), Rk,e(t) indicate the densities 
of S-, I-, R-nodes at time t, respectively. Naturally, 
by setting k = £, the same arguments are applied to 
undirected networks. 



A. Imperfect immunization to targeted nodes 

This is the case where all of target nodes in T cannot 
be immunized, because some are overlooked or hidden 
etc., and a is like the vaccine coverage within T. 

Since there are no rewiring of links and immunized 
nodes are no longer S-nodes, we have 



SkAt) + PkAt) + RkAt) 



l-a, (k,£)GT, 

1, (k,£)&T, W 



The condition a = 1 implies the perfect immunization 
for T, while a = means no immunization, which co- 
incides with the simple SIR model [11]. 

Following the dynamical mean- field theory ([6, 9, 
10]), we see that the spreading process on a directed 
network can be described by the system of differential 
equations: 



dSki 

dt 
dpk,£ 

dt 
dRkj 

dt 



-\kS k At)0{t), (2) 
XkS k At)0(t) - p k At), (3) 
Pfc,/(*)- (4) 



The term XkSk,e(t)8(t) in (2) and (3) indicates the frac- 
tion of newly infected nodes through fc incoming links, 
while 9(t) is the probability of contact with I-nodes 
from whose outgoing links the disease spreads. Hence 
6{t) can be written in parts as 



+ J2£p(k,£)p k At)), (5) 



for the later use. Here the denominator is the average 
{£) and each directed link is counted twice as one out- 
degree of some node and as one indegree of another. 
Hence the average outdegree is equal to the average 
indegree: {£) = (k). From (1) the initial conditions 
imposed upon Sk,e are 



S fc) /(0) 



l-a, {k,£)eT, 
1, (k,£)ef. 



(6) 



The initial values of other variables are assumed to be 
zero or nearly zero. 

B. Insufficient immunity 

This model is applied to the case, for example, 
where vaccination or a treatment cannot provide suf- 
ficient immunity upon the target nodes in T, so the 
rate a (0 < a < 1) represents the immunity effect. The 
higher it is, the less infected S-nodes in T become. The 
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differential equations are divided into two parts. First, 
for (k,£) e T, it follows that 



dSkj 
dt 

dt 
dRkj 
dt 



-Xk(l-a)S k , t {t)0(t), (7) 
Afc(l - a)S k j{t)6(t) - p k>e (t), (8) 

PkA*)- (9) 



Secondly, for (k, I) G T, Eqs. (2)-(4) remain valid. The 
probability 9 is the same as (5) and the initial condition 
imposed upon S k ,e is 

ShAo) = i 

for all k and I. The initial values of other variables are 
zero or nearly zero. 

The next section analyzes models A and B in di- 
rected complex networks with target set T and with 
immunization rate a. In Section 4 the corresponding 
SIS models are considered. Actually all thresholds turn 
out the same. It is easy to see that the following argu- 
ments can be applied to models A and B considered in 
undirected complex networks. 



3. Epidemic spreading in the SIR model 



A. Imperfect immunization to targeted nodes 

To begin with we solve Eq. (2) under the initial 
conditions (6): 

f (l-a)e-***<*> for(M)€ET, 
S MW- e -Xk4>(t) t^rue\^m ,l! 



for (7c, f) € T, 



where 



<f,(t) = f 6(t')dt'. 
Jo 

By (4), (5) and R k ,e{°) = for a11 k and ^> has an 
expression: 

+ f PkAt'W) (n) 

f Jo 

+ ^ P (M)fl M (t)}. 

T 

To derive the differential equation for <p(t), we use (1), 



(10) and (11). Then it follows that 
d(j){t) _ 1 



dt (£) 



— { e m*. + E ^( fc ' } 



i-{EMM)(i - " - fl M (*) - s M (t)) 



A]^,Q(i-fi M (f)-S w (f))} 



= ^{£MM)(i -<*-&,/(*)) 

+ -&,/(*))}-#*) (12) 

T 

= £MM)(l """(I- «)e- Afc0(t) ) 
f 

= ^{(i-a)£W)(i-e-^«) 

+ EM*,*)(l-e- Afc * (t) )}-0(t). 
f 

We are concerned with a steady state of the epi- 
demic spreading, for which one has a limit 

<f> = lim (j){t), 

t— >00 

together with the condition 

lim«=0. 

t >oq at 

Substituting these two into Eq. (12), we get the equa- 
tion for $ as follows: 



$=i{(l- a )^#,f)(l- £ - M ) 



(13) 



An epidemic outbreak implies that this equation has 
a solution $ > other than $ = 0. Since the right hand 
side of (13) is a monotone increasing and concave func- 
tion of $ and its value at $ = 1 is less than 1, the 
condition for it is 

l^{(l-")EW)(l-e-^) 

+ EMM)(i- e - Afe *)}L =0 >i. 

T 

Thus the critical threshold A c is expressed as 



A c = 



(1 - a) J2t ke P( k > l ) + Et k£ P(k, i) ' 
Alternatively, using the identity 

E^P(M) = (k£) -^2k£p(k,£), 



(14) 
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(14) is equal to 



we can write the fraction (16) of the outbreak size as 



(k£)-a(k£) T 



(15) 



In [7, Chap. 10] the total number of infected indi- 
viduals is discussed for the classical SIR model, which 
represents the final outbreak size O. In our setting it 
is the averaged fraction of nodes ever infected until the 
disease dies out. Using the solutions (10), this is writ- 
ten as 

= ]>>(M)(i-SmM) 
k.t 

= 1 - 5>(M)S fcl /(oo) - $>(M) S m(°°) 

T f 

= 1 - (1 - a) ]>>(MK Afe * - 5>(M)e" Afc * 

T f 

= i - £>(M)e- Afe * + «]>>(M) e - Afc$ 

= l-^P(fc) e - Afe *+ a ^p(A : ,^) e - A ' £ *. (16) 



The last term of the right hand side provides the effect 
of an immunization. 

As an illustrative example let us consider an im- 
munization of nodes with large indegrees. This means 
T={(k,t)\k > K} or 

Sk,t(t) + pkA t ) + R kA t ) = | \ a ' \ < k' 

where if is a certain large value of the indegree. The 
corresponding initial conditions imposed upon Sk,e are 

. . / 1 - a, k> K, 
Sk ^°> = { 1, k < K. 

The initial values of other variables are zero or nealy 
zero. 

Using the solution $ of Eq. (13) for this case, we 
obtain the fraction (16) of the outbreak size: 



0=l-J2 P{k)e- Xk * +aJ2 p ( k ) e 



\k<S> 



k>K 



by means of the indegree distribution P(k). If the in- 
degree distribution follows a power law 

P(k) oc k ~ 7 , k > to 

with 7 > 1, then P(k) = (7 — l)m 7_1 fc -7 (k > to), 
where to is the minimum indegree. Defining the incom- 
plete gamma function [1, p. 260] by 



POO 

T(a,x) = / f-V**, 

J X 



O 



P OO f-CO 

= 1- P(k)e~ Xk *dk + a / P(k) 

Jm J K 



e xk *dk 



= 1 - (7- l)(Am$) 7 - 1 r(l - 7 ,Ato$) 
+ a(7 - 1)(XK $) 7 " 1 r(l - 7, \K§), 

by the continuous approximation. 

B. Insufficient immunity 

Since similar arguments used in A are also applied 
to this model, we only give an outline of the analysis. 
Solving (7)-(9) under the initial condition 5^,^(0) = 1 
for all (k,£), we have, in place of (10) and (13), 



Sk,e{t) 



e -X(l-a)k<P(t) { 0T (k,e)€T, 

e -\k4>(t) {or(k,l)€T, 



and 



$ = i{^MM)(i-^ AMB ) 



+ ^Mfc^)(l-e- Afc *)}, 



respectively. Solving the following equation for A, 
^{£MM(l-e- A(1 - Q)fc *) 

+ 5>(M)(i- e - Afc *)}| $=0 = i, 

T 

we again obtain (14) and (15). When T = {{k,t)\k > 
K} as above, the averaged fraction of nodes ever in- 
fected until the disease dies out becomes 

O = 1 - ^ P(k)e- x(1 - a ^ - P{k)e- Xk *, 

k>K k<K 

by an application of (16). 

4. Epidemic spreading in the SIS model 

In the SIS model, R-nodes are absent and nodes 
that recovered from the disease cannot acquire eternal 
immunity and may be infected again and again. Let T 
and a (0 < a < 1) be the same as before in the mod- 
els A and B. The differential equations for the densities 
Sk,i(t) 7 Pk,e(t) of S-, I-nodes at time t are easily derived 
from those for the SIR models in Section 2. 

A. Imperfect immunization to targeted nodes 

In this case, remarking (1), we have 

f l-a-p M (i), for (h,t)€T, 

bkA j I 1-PkAt), for(M)eT. 
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Therefore, Eqs. (2)-(4) in Section 2 are replaced by the 
differential equations 



dt 



' Afc(l-a-pfc,/(t))0(t)-p fc ,/(t) ) 

if (k,£) gT, 

if {k,£) e? 1 , 



(17) 



where is the same probability as (5). 

At the steady state, as in Section 3, we will have 
the condition 



lim 

t-s-oo at 

for all k and £, and a limit 







6 = lim 9{t). 

t—t>oo 

So we get from (17), 

f (l-a)Afc©/(l + Afc©), if (M)eT, 

mnp MW -| Afce/(1 + Afce); if(fc,*)eT. 

Substituting these into (5) and letting i — > oo, we 
have the equation for as follows: 



e = 



^A/Wp(M)©)- 



+ 



If this has a solution > other than = 0, then it 
corresponds to an endemic state. Since the right hand 
side of the equation is a monotone increasing and con- 
cave function of and its value at = 1 is less than 
1, the condition for an endemic outbreak is 

d 1 t 1 ~ u)\k£p{k,£)Q 



dQ (£) V 



^ Xk£p(k,£)Q 
\* 1 + Afc0 



1 + A/J0 

) 



> 1. 



e=o 



Again, this yields the same threshold (14) or (15) as 
the model A for the SIR model: 



{k£} -aJ2 T ki P( k ^) (k£)-a(k£) T 
B. Insufficient immunity 



This model just coincides with Eq. (10) of [4] for 
the SIS model on undirected networks. Taking Eq. (3) 
for T and Eq. (8) for T, we get 



dt 



{ Xk(l - a)(l - p k , e (t))6(t) - p k ,t(t), 

if (k,£) gT, 

Xk(l- Pk At))0(t) - p k At), 

if (k,£) GT. 



Using this equation for (17) and repeating the above 
calculations, we obtain again the same threshold (14) 
or (15) as in A. 

5. Immunization strategies 

We have seen that the thresholds are all the same 
in four cases. So in any one of them we can evaluate 
immunization strategies based on the thresholds. In 
view of (14) or (15), to control epidemic outbreaks, it 
is desirable to choose the set T so that (k£) T becomes 
as large as possible, and T should be a considerably 
small region compared to the entire population from a 
practical point of view. 

First remark that in case of no immunization (a = 
0) the epidemic threshold is equal to 



A° = 



(kty 



by (14) or (15). 

If the entire population is an object for immu- 
nization, that is, T is empty, then from the identity 
(k£) T = (k£) the threshold due to this global immu- 
nization is 



A r 



(i) 



1 



(k£) - a{k£) T 1 - a (k£) ' 



Suppose a < 1 . Then the threshold for any targeted 
immunization satisfies 



A c = 



(k£) - a{k£) T 



< A c = 



1 



1 - a 



(18) 



since (k£)r < (k£) always holds. Therefore, in case of 
a < 1, certainly any targeted immunization is not more 
than the global immunization. 

Suppose that in addition to a < 1 considered net- 
works are scale-free as presented in Section 1. Scale- 
free networks are prone to have very large (k£); A° w 0, 
under a high correlation between nodal indegrees and 
outdegrees [11]. In such a case inequality (18) tells us 
that any targeted immunization is not effective: A c w 0. 

Next supposing a = 1, let us consider a simple tar- 
get set T by defining T = {(k,£)\k,£ < M} for some 
positive M. This yields 

(k£) f = k£p(k,£) < M 2 , 

kJ<M 

and hence, even in scale-free networks, we are able to 
conclude 



(k£) - (k£) T (k£) 7 



{(■) 



M 2 



> > K « o, 



which implies the effectiveness of the targeted immu- 
nization only under the condition a = 1. The choice 
of a desirable target set T depends on the probability 
distribution p(k,£), so this estimate may be improved 
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further. 

Finally, as for undirected networks, there is no dis- 
tinction between indegrees and outdegrees. So, setting 
k = £ in (15), the threshold of an epidemic outbreak 
can be written in our notation as 

c (fc 2 ) -a(jfc 2 ) T ' 

for an immunization rate a and a target set T that is 
characterized by the degree k. Using this, we see that 
similar results of this section hold for undirected net- 
works. Note that in [12] the target set was assumed to 
be T — {k\k > K} for a large number K, namely a set 
of hub nodes. In this case the targeted immunization 
is not effective unless it confers 100 percent immunity 
on all hub nodes in T. 
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